Can parameters of / - mesons t> e determined correctly analyzing only tytv scattering? 
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The coupled processes - the mv scattering and tvtt — > KK in the I G J PC = + ++ channel - are 
analyzed (both separately and combined) in a model-independent approach based on analyticity and 
unitarity and using a uniformization procedure. It is shown: 1) a structure of the Riemann surface of 
the S'-matrix for considered coupled processes must be allowed for calculating both amplitudes and 
resonance parameters, such as the mass and width; 2) the combined analysis of coupled processes 
is needed as the analysis of only mr channel does not give correct values of resonance parameters 
even if the Riemann surface structure is included. 
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\ The study of scalar mesons is very important for such profound topics in particle physics as the QCD vacuum. 
However, despite of a big effort devoted to studying various aspects of the problem [l[ (for recent reviews see 
a description of this sector is far from being complete. Parameters of the scalar mesons, their nature and status of 
' r~ I '■ some of them are still not well settled m. E.g., applying our model-independent method in three-channel analyses 
O I of processes nir —> 7T7T, KK, T]T](T]r]') @,T3] we have obtained parameters of the /o(600) and /o(1500) which differ 
' . considerably from results of analyses utilizing other methods (mainly those based on the dispersion relations and 
Breit-Wigner approaches). Reasons of this difference should be understood because our method is based only on 
! a demand for analyticity and unitarity of the amplitudes using a uniformization procedure. The construction of 
amplitudes is practically free from any dynamical (model) assumptions using only the mathematical fact that a local 
behavior of analytic functions determined on the Riemann surface is governed by the nearest singularities on all 
corresponding sheets. I.e., the obtained parameters of resonances can be considered as free from theoretical prejudice, 
^v^j ' First note that in our previous three-channel analyses with the uniformizing variables @, @] we were enforced to 
(X) ■ construct a four-shcctcd model of the cight-shccted Riemann surface. This we have achieved neglecting the tttt- 
lO ' threshold branch-point which means that we have considered the nearest to the physical region semi-sheets of the 
initial Riemann surface. This is in the line with our approach of a consistent account of the nearest singularities 
on all relevant sheets. The two-channel analysis utilizes the full Riemann surface and is, therefore, free of these 
approximations. To verify a plausibility of our assumptions in the three-channel calculations, we have performed a 
■ combined two-channel analysis of data on irir — > ttit, KK to check whether the results of our three-channel analyses 
[1, 0] are also obtained in the two-channel consideration. Moreover, to better understand reasons for the above- 
indicated difference in results, we have performed first the analysis only of the 7T7t scattering data in the two-channel 
approach. 

rS \ The two-channel S'-matrix is determined on the four-sheeted Riemann surface. The matrix elements Sij, where 
i,j = l(mr),2(KK) denote channels, have the right-hand cuts along the real axis of the s complex plane (,s is the 
invariant total energy squared) , starting with the channel thresholds Si , and the left-hand cuts related to the crossed 
channels. The Riemann-surface sheets are numbered according to the signs of analytic continuations of the roots 

^Js — Si as follows: signs (im^/s — s\, Im-^/s — 82) = ++, — +, , H — correspond to sheets L II, III, IV. 

The resonance representations on the Riemann surface are obtained from formulas |8| which express analytic 
continuations of the S'-matrix elements to unphysical sheets in terms of those on the physical (I) sheet having only 
resonance zeros (beyond the real axis). Then, starting from the resonance zeros on sheet I, one can obtain an 
arrangement of poles and zeros of a resonance on the whole Riemann surface ( "pole clusters" ) . In the two-channel 
case, according to these formulas and a real analyticity amplitudes of all coupled processes have conjugate poles at 
the same points of complex energy on sheets II, III and IV if S^ = 0, S^S^ — (S{ 2 ) 2 — an d S\ 2 = 0, respectively 
(the superscript I means the matrix elements on sheet I). Therefore three types of pole clusters representing states 
of different nature arise: (a) when there is a pair of conjugate zeros only in S\ x = 0, (b) only in S\ 2 = 0, and (c) in 
both the S\ x = and S\ 2 = 0. If the coupling of channels is present (S12 ^ 0) (i.e. a state decays into both channels 
and/or is exchanged in crossing channels), then this state being of type a is represented in Sn by a pair of conjugate 
poles on sheet II and a pair of conjugate zeros on sheet I and also by a pair of conjugate poles on sheet III and a pair 
of conjugate zeros on sheet IV at the same complex-energy points, which are to be shifted with respect to the zeros on 
sheet I. For states of type b, the pair of conjugate poles on sheet III is shifted relative to the pair of poles on sheet IV. 
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For the states of type c, one must consider two pairs of conjugate poles on sheet III. Generally, wide multi-channel 
states are most adequately represented by pole clusters, because the pole positions are rather stable characteristics for 
various models, whereas masses and widths are very model-dependent For calculation of the latters one must use 
the poles on those sheets where they are not shifted (due to the channel couplings) with respect to the zero position 
on sheet I. For resonances of types a and b these poles are on sheets II and IV, respectively. For resonance of type c 
the poles can be used on both these sheets. In the case of N channels, the poles only on the sheets with the numbers 
2* (i = 1, ■ ■ ■ , N is the number of channel), i.e. II, IV, VIII,. . ., should be used for calculating resonance parameters 
0,|. 

It is convenient to use the Le Couteur-Newton relations [K|. They express the S- matrix elements of all coupled 
processes in terms of the Jost matrix determinant d(yj s — s\, • • • , \/ s — sjv) that is a real analytic function with the 
only square-root branch-points at y/s — s a = 0. 

A necessary and sufficient condition for existence of the multi-channel resonance is its representation by one of the 
types of pole clusters. To use this representation of resonances, which is very important for the wide multi-channel 
states, a uniformizing variable is applied. Analyzing hit — > nn, KK we applied the uniformizing variable Q which 
takes into account, in addition to the irir- and ATA'-threshold branch-points, the left-hand branch-point at s = 0, 
related to the irn crossed channels: 



rriK s/s — Am% + m w y/ s — 4m 
V s ( m K ~ m l) 



1 (1) 



It maps the four-sheeted Riemann surface with two unitary cuts and the left-hand cut onto the w-plane. Representation 
of resonances of various types on the uniformization w-plane can be found in [ill ] . 

On the v-plane, Sn(v) has no cuts; S 2 2 {v) & n d S22{v) do have the cuts which arise from the left-hand cut on the 
s-planc, starting at s — A(m 2 K — to 2 .), which is further approximated by a pole 

d L = v- 4 (l - (p - i^/l - p 2 )vY(l + (p + i^l- p 2 ) 



where p = 0.903 ± 0.0004 from analysis. An explanation of the fourth power of this pole can be found in [Tl| . 

On the f-plane, the function d(v) in the Le Couteur-Newton relations [1, [To| does not possess any branch points. 
The main model-independent contribution of resonances, given by the pole clusters, is factorized in the S'-matrix 
elements from the background. The possible remaining small (model-dependent) contributions of resonances are 
supposed to be included in the background. Therefore, d(v) is taken as d — d res dLdb g where the resonance part is 



d res = v~ M Y[(l - v*v)(l + v n v) (2) 



n=l 

with M the number of pairs of the conjugate zeros. The background part is 
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d bg = expj^-zy^(Vs - s n /2m n )(a n + i[3 n ) , (3) 



n=l 



a n ,0 n = a n i,b n i + 2^ a n k,b n k(s/sk — l)0(s — Sk) 

k=rj 7 a,v 

with s v and s a the rjrj and era thresholds, respectively, s„ a combined threshold of the 7777', pp and ujuj thresholds; 
from the analysis: s a = 1.6558 GcV 2 , s v = 2.1293 GeV 2 . 

At present in the scalar sector, there are alternative data for the nir scattering - [l2[ and [HI - which arc different 
considerably in the 0.76 GeV region and especially above 1.45 GeV. Therefore separate analyses using these alternative 
data are needed. Here we performed the analysis taking for the irir scattering in interval 0.575 GeV< v / s< 1-89 GeV 
the data from [l2| and for ^fs < 1 GeV from many works. References to the latter and to practically all accessible 
tttt — > KK data used can be found in [Tl| - Analysis using the data [l]| will be presented in other paper. 

First analyzing only the tttt scattering, we supposed an existence of two states (narrow /o(1500) and wide /q(1500)) 
and achieved an excellent description for the phase shift Sn and modulus |5n| (the total x 2 /NDF« 1.07) with the 
resonance parameters (Table HJ which largely coincide with estimations of the PDG [l| (cf. also [HI)- The only 
distinction is observation of the wide /q(1500). The fact that this state is not observed in works cited by the PDG is 
related, it seems, with peculiarities of analysis of data therein. 
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TABLE I: The pole positions of resonances ^/s^ = E r — iF r /2 [MeV] on sheets II and IV in the y^s-plane in the analysis of only 
7T7r-scattering. 



Sheet II IV 

/o(600) 447.5±5.9-i(267±6.5) 
/o(980) 1001.1±3.7-i(20.3±2.6) 

/o(1370) 1301.1±47.9-i(224±49.3) 
/o(1500) 1503.7±45.1-i(56.5±39.4) 
/o(1500) 1511.4±11.2-i(200.5±ll) 1505.9±38.5-i(168±40.6) 
/o(1710) 1720±32.2-i(64.9±30.1) 




FIG. 1: The S-wave phase shifts and modules of the 7T7r-scattering and tttt — > KK matrix elements. The short-dashed lines 
correspond to the analysis only of the tttt scattering. The long-dashed and solid lines correspond to solutions A and B of the 
combined analysis of tttt tttv, KK, respectively. The data are from Refs. [HI . [T3 , fTH. [TB|. 



In the analysis, the /o(600) and /o(980) are described by the clusters of type a; /o(1370), /o(1500) and /o(1710), type 
b; /^(1500), type c. The received background parameters are: on = -0.0895±0.0030, ai, = 0.04±0.03, ai CT = 0.0±0.8, 
ai^=0.0±0.7, &u=0.0±0.007, b lv = 0.0±0.01, b la = 0.0±0.02, bi„=0.054±0.036. 

In Fig.Q] we show the fitting only to the tttt scattering data and energy behavior of the phase shift 0i2 and modulus 
of 5*i2 calculated using the resonance parameters from this analysis. In spite of the very good description of data 
and the very good agreement of obtained 7T7T scattering length a® with the experimental results and with the chiral 
perturbation theory (ChPT) calculations (see Table iHTj) . this analysis shows two important flaws: (1) The negative 
background phase-shift beginning at the tttt threshold (an = —0.0895) is necessary for a successful description of the 
data. This should not be the case because, in the uniformizing variable, we have allowed for the left-hand branch-point 
at s = which gives a main contribution to the tttt background below 1 GeV. Other possible contributions of the 
left-hand cut from exchanges by the nearest p and Jo (600) mesons practically obliterate each other Q because vector 
and scalar particles contribute with the opposite signs due to gauge invariance. (2) Description of the 7T7t —> KK 
data, using the same parameters of resonances as in the 7T7t channel, is satisfactory only for the phase shift <pi2 which is 
due to the approximation of the left-hand cut in S\2 and £22 by the fourth-power pole. The modulus |5i2| is described 
well only from the KK threshold to about 1.15 GeV as it should be due to the two-channel unitarity. Above this 
energy the description fails even qualitatively (Fig. [I}. 

From this we conclude: If the data are consistent, for obtaining correct parameters of wide resonances the combined 
analysis of coupled processes is needed. Further that analysis of 7T7t — »• 7T7t, KK is performed successfully. The data 
for the 7T7T scattering below 1 GeV admit two solutions for the phase shift - A and B - which differ mainly in the 
pole position on sheet II of the /o(600). The total x 2 /NDF is 1.53 for the A-solution and 1.44 for B-solution. The 
resonances are described by pole clusters of the same types as in the analysis only of the tttt scattering. In Table |TT] 
we show the pole positions of resonances on sheets II and IV on the \/s-plane. 
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TABLE II: The pole positions of resonances on sheets II and IV in the y^i-plane in the combined analysis of the -kit —¥ tttv, KK 
data. The complete pole-clusters of resonances can be found in [llj . 



Sheet 


II IV 


A-solution 


/o(600) 


517±7.8-i(393.9±6) 


/o(980) 


1004.6±3.9-i(25.0±2.3) 


/o(1370) 


1342.9±12.2-i(221.6±30.7) 


/o(1500) 


1501.1±6.4-i(56.6±6.0) 


/o(1500) 1532.2±12.4-i(323.2±21) 1519.3±18.7-i(339.5±42.2) 


/o(1710) 


1717±34.9-i(72.9±16.2) 


B-solution 


/o(600) 


550.6±9-i(502.1±7.2) 


/o(980) 


1003.2±3-i(28.9±2) 


/o(1370) 


1336.7±14-i(251.9±27.5) 


/o(1500) 


1500.3±6.3-i(57.0±6.4) 



/o(1500) 1528.4±12.5-i(328±20.2) 1515.6±17-i(340.3±34.9) 
/o(1710) 1722±35.7-i(92.3±20.3) 



The obtained background parameters for the A-solution arc: an = 0.0 ±0.003, cii n = — 0.1004 ± 0.0301, ai a = 
0.2148±0.0822, a lv = 0.0±0.07, b n = b lv = b la = 0, b lv = 0.012±0.0287, a 2 i = -0.919±0.107, a 2n = -1.399±0.348, 
a 2cr = 0.0±0.7, a 2 „ = -11.45±0.75, b 21 =0.0747± 0.05 03 , 6 2t? = o 2(T = 0, b 2v =4.83±1.94; for B-solution: a n =0.0±0.003, 
ai, ; = -0.0913±0.0327, ai CT = 0.1707±0.0899, a lv =0.0±0.07, b n = bir, = b lcr = 0, b lv = 0.006±0.029, a 21 = -1.338±0.111, 
a 2t; = -1.119±0.376, a 2a = 0.0±0.8, a 2v = -12.13±0.77, b 21 =0.018±0.050, 6 2r; = fe 2(T = 0, b 2v =4.48±1.98. 

In the combined analysis both flaws of the only 7T7r-scattcring analysis arc cured. Now the 7T7t background below 
the KK threshold is absent (an = 0.0). An arising pseudo-background at the r\r\ threshold (ai,, < 0) is also clear: 
this is a direct indication to consider explicitly the ry^-threshold branch-point. This was already done in our work 
0. In the combined analysis the /o(600) parameters are changed considerably receiving new values closer to those 
obtained in our three-channel analysis [7J. Earlier one noted that wide resonance parameters are largely controlled 
by the non-resonant background [171 ] . In part this problem is removed due to allowing for the left-hand branch-point 
at s = in the uniformizing variable. 

In the Table IIIII we compare our results for the 7T7t scattering length a[J with results of some other theoretical and 
experimental works. 



TABLE III: The tttt scattering length ag. 



a° [m-l] 


Remarks 


References 


0.222 ± 0.008 


Analysis only of 


This work 




7T7T scattering 




0.230 ± 0.004 


A-solution 


This work 


0.282 ± 0.003 


B-solution 


This work 


0.26 ±0.05 


Analysis of the K — > -nirev 


[15J 




using Roy's equation 




0.24 ±0.09 


Analysis of -n~p — > 7r + 7r~?i 


[16j 


0.2220 ± 0.0128 stat 


Experiment on K e 4 decay 


[18] 


±0.0050 syst ±0.0037 th 






0.220 ± 0.005 


ChPT ± Roy's equations 


[19, 20] 


0.220 ± 0.008 


Dispersion relations 


[21] 




and K e 4 data 




0.26 


NJL model (I) 


[22] 


0.28 


NJL model (II) 
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In the analysis only of tttt scattering and in the A-solution we reproduced with a high accuracy the ChPT results 
fl9l l20j including constraints imposed by the Roy equations. On the other hand, the B-solution is similar to the 
predictions of the chiral approaches based on the linear realization of chiral symmetry (models of the Nambu-Jona- 
Lasinio (NJL) type [22|, [23 )■ Taking into account very precise experiments at CERN performed by the NA48/2 (l8j 
and the DIRAC [24[ Collaborations, which confirmed the ChPT prediction [l9|,[20|], one ought to prefer the A-solution. 

In summary, a structure of the Riemann surface of the 5-matrix for coupled processes must be included properly. 
To calculate resonance parameters, such as masses and widths, one must use poles on those sheets where they are not 
shifted (due to the channel coupling) in respect of the zeros on sheet I. In the two-channel case, the relevant poles 
are on sheets II or/and IV depending on the resonance type. Moreover, the combined analysis of coupled processes 
is needed as the analysis of only 7T7t channel do not give correct values of resonance parameters even if the Riemann 
surface structure is included. Finally, in order to be concrete, the main scope of the paper is the scalar mesons, 
however, the method and conclusions of paper can be applied to other wide resonances, e.g. to vector mesons. 
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with JINR, the Polish Ministry of Science and Higher Education (grant No N N202 101 368) and by the DFG under 
Contract No. LY 114/2-1. 
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